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ABSTRACT 


An opacity sampling (OS) technique for treating the radiative 
opacity of large numbers of atomic and molecul.tr lines In cool stellar 
atmospheres Is subjected to several tests. In accord with previous 
results, we find that the structure of armosphe *lc models Is accurately 
fixed by the use of 1000 frequency points, and S00 frequency points Is 
often adequate. The effects of atomic and molecular lines are 
separately stidled, A test model computed using the OS method agrees 
very well with a model having Identical atmospheric parameters computed 
with the giant line (opacLty distribution function! method. 
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I. INTRODUCTION 

An opacity sampling (OS) technique appears promising for treating 
the opacity of large numbers of overlapping molecular and atomic linen In 
the spectrum of coal stars. In the OS method for model atmospheres, no 
approximation whatever Is nadj In the calculation of the atomic and molecular 
opacity at a selected frequency. Hie resulting flux (F(v, depth)) consequently 
contains no errors due to opacity smearing. In a computation of a model 
atmosphere, then, the only approximation In computing the opacity Is In the 
choice of a particular finite set of frequency points. (The atmosphere 
would still he subject to such other approximations as plane parallel geometry 
and LTE, of course.) 

The method Itself has been derived Independently and applied to the 
treatment of atomic lines In the Sun (Peytremann 1974) and to molecular lines 
In the atmosphere of a carbon star model (Johnson and Sneden, 1974; Sneden, 
Johnson, and Krupp, 1975 hereafter referred to as Pnpev I). These earlier 
papers have discussed the advantages and dl sadvantages of th“ OS method 
compared to other line blanketing methods. 

In the present paper we apply new, stringent test9 to the self-consistency 
of the method and determine the sensitivity of the resulting model atmospheric 
6truct'jre to the number of frequency points used. We then extend the treatment 
to Include the combined effects of atomic and molecular line blanketing. We 
apply the method to the calculation or model atmospheres for a K giant star and 
compare our results to those obtained using either straight mean opacities or 
opacity distribution functions (ODF) . 


II. TECHNIQUES ANT) DATA 


Since the basic philosophy and a quantitative description of the method 
of opacity aampllnR has already been given In the literature (Peytremnnn 1974; 

Paper I), w lve here only a brief summary of the 

aLthod. The method Itself Is th». sence of simplicity. At a selected set 
of frequencies (v) and optical depths (t) or maos points (m) one calculates 
the correct linear or mass absorption coefficient by adding together all the 
contrlbr o..s of all the atomic and molecular lines as well as the continuous 
sources which contribute at that wavelength. The number of such lines may 
range from, a few to a few hundred at each frequency. Such a method obviously 
requires an enormous data bank, but these are now becoming available. Once one 
has an accurate value for the absorption coefficient K(v,m\ one can calculate 
every monochromatic quantity Q(v,m) needed to produce the model atmosphere. 
Monochromatic fluxes ctmputed bv this method should be perfectly accurate, 
within the usual framework of plane-parallel geometry and LTE. Energy fluxes 
Integrated over the entire spectrum should be exact, the onl" app-oxlmat ion being 
In the number of frequency pol'ts. Calculations have already Indicated that for 
either atoms or molecules the rnumber of frequency points needed Is between 500 
and 1000 (Peytremanu 1974; Paper I). We Investigate 

thin question further In the present paper. 

To compute model atmospheres we use a version of the compute; code ATLAS5 
(Kurucz 1970), which we have modified to Include monochromatic molecular and 
atomic line absorption. The code Is based on hydrostatic equilibrium and constancy of 
total energy flux (radiative plus convective) In a plane-parallel atmosphere. 

We make the usual assumption thee the energy level populations of both atoms and 
molecules are distributed according to LTE and assume that all lines are formed 


in pure .'b'nrptlon. 
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The line data was assembled from several sources. As part of a long 
term project, we have previously assembled at Indiana University several data 
tapes containing Information on over one million lines of Cp (^wan, Phillips, 
and Balllk-Ram-iey systems), CN (red and violet systems), CO (Infrared), and 
CH (G-band system). The wavelengths and line strength factors for these lines 
were obtained from computer programs written by T, D. Fay and based on the 
formalism of KovacB (1969) , using the best available (-values and constants in 
the literature. Roger Bell kindly provided us with a data tape containing 
empirical Information on approximately 35,000 atomic lines as well as several 
thousand lines of MgH, NTH, and OH. Finally, In the Intervals X < 3000A and 
\ > 6500A we Include Fe I and Fe II lines from u semi empirical data tape 
provided by Rober, Kurucz (1974). 

For both the tests of internal consistency and comparisons with other line 
blanketing methods, only the Doppler core of the line absorption coefficient 
Is considered. The addition of natural and collisions! damping to the line blanketing 
has -'egligible effect In the models used here. The Ca II H and K lines 

are not included In these calculations because their great strength would 
necessitate special treatment, and, for these cool stars, only a small fraction 
of the energy flux is contained In tne spectral interval X * 4000A, 

To prepare the line data for use In the model atmosphere program we first 
select a set of frequencies at which the radiation field la to be calculated. For each 
of these frequencies, we abstrac from the atomic and molecular tapes 
all those lines which fall within a given Interval around the selected frequency. 

We commonly use the Interval Av ■ ± 1cm * for molecular lines and twice that 
value for atomic lines. All lines are then sorted 


by wavenumber and stored on a tape ready for use In the mod* 1 atmosphere program. 



The various test models were normally converged until a relative flux 
constancy to within 0.0005 was achieved; at this point the Indicated temperature 
corrections were usually smaller than 2*K. For a model started from an 
ordinary guess, approximately 20 lreratlons were required to restore flux 
constancy, requiring 1 hour on a CDC 6600 If 1000 frequency points were used. 
Faster convergence can be achieved by occasional "hand" smoothing of the 
Intermediate model atmospheres, as well nt by choosing a better "guess" as to 
the final model. 


III. RESULTS AND DISCUSSION 


A. Internal Consistency 

For the same set of physical conditions we compute models with differing 
numbers of frequency points and spacing* In order to determine the number of frequencies 
are req .ired to compute accurate model atmospheres self-consistent In such a way 
that: 1) a further Increase In the number of frequencies does not appreciably 

change the resulting atmospheric structure, and 2) for that critical number of 
frequencies, the atmospheric sfucture Is relatively Insensitive to the specific 
frequency points chosen. 

Our standard model has the following characterise ^ ff " K; 

2 

surface gravity ■ 100 cm/sec'; v (uiic ro turbulence) ■ 1 km/ sec (used only in the 
line absorption coefficient, not the pressure); solar abundances (Lambert 1968); 
continuous opacity sources as In ATLAS 5; and line opacities treated by the OS 
method. This model will be referred to as (4000/2. 0& /HOP, OS). 

For purposes of Illustration a standard modal a'-iiMsphere using 1010 
frequency points has been adopted (which we, for our convenience, will always call 
a 1000- frequency set). This number Is suggested by previous studies in which 
500-1000 frequencies (Taper I) or 600 (Peytrem.ann, 1974) frequencies wen found 
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adequate to apeclfy accurately th.- structure of a model atmosphere. The 1000 
frequency set Is equally spaced by 25cm 1 from 500 to 25,000 cm* 1 (20 u to 

• . j 

4,000a); an equal spacing of 250 cm la used \n the ultraviolet down to a 
wavelength of 27004. A coarser grid of points Is employed In the region 
2700 < X < 4,000A because this portion of the spectrum transmits only a small 
portion of the total flux In our models, and the flux shortward of 2700A Is 
totally negligible for photospherlc models at these low temperatures. 

To test the effects of the frequency set on the model atmosphere, the 
following grid of standard models has been calculated: three models with different 

sets of 100 frequencies, three models with different sets of 500 frequencies, 
two models with different sets of 1000 frequencies, and one model with set of 
2000 frequencies. 

The three 100-frequency sets (RlOlA, F101B, R101C) are semi -random in their 
spacing, obtained from a 1000 point set similar to the one described above. 

That la, one point in every ten Is chosen from the larger set. with the help 
of a random number table. This procedure guarantees 

an appropriate dispersion of points across the entire spectrum, yet each frequency 
set so generated Is completely Independent of the others. 

Two of the 500-frequency sets (500A, 500B) are spaced equally by 50cm * 
below 25,000 cm , ana above by 500 cm *. The two sets are shifted with respect 
to each other by 25cm *. The third (R500) ».s made up of random frequencies; 

Co Insure coverage of the entire spectrum, however, one point Is randomly selected 
from each Interval of 100cm * for wavenumbers below 25,000cm * (4,000A) and In 
each Interval of 1000 cuT* for wavenumbers above 25,000 cm'*. All three sets are 
equally unbiased with respect to the positions of observed spectral lines. 

One of the 1000 frequency sets Is the standard set described above: the second 

-1 -1 

has a frequencies randomly selected within 50cm bins (500cm bins for wave- 
numbers i.liove 25,000 cm *) across the spectrum. 
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The result* of these model atmosphere calculations are shown In Figures 
1 and 2, In which .ire plotted the temperature versus th* mass parameter 

r\ , p 

RHOX ■ iqdx, which represents the mass In a unit column (g cm ) above that point. 

Since total pressure ■ surface gravity x RHDV, these graphs are T-P telatlon* 
at the same time. Figure 1 shows the three 100- frequenc y point models plotted 
together with the 1000-polnt standard model, represented by a solid line. 

Vertical bars Indicate a few selected values o f log r (where t, ■ Rosselnnd optical 
depth). At RUOX • 1 there la a spread of 2S0*K among the various models, showing 
that the upper atmosphere It- no' well represented by a 100-polnt frequency set. 

The spread Is even greater at smaller values of RHOy, with a maximum deviation of 
about iy # K. However, the photospherlc ' emperaru re -pressure structure (in the 
region 0.01 4 t £ 5, where the continuum Is formed) Is matched quite well. 

ft 

Figure 2 Is a similar diagram Involving the three *>00 frequency point 
models. At RHOX * 1 there Is a spread of approximately 100*K among the various 
models, a value substantially smaller than In the previous case. A temperature 
difference of this order would he marginally detectable by observations. Once 
again, the photospherlc structure Is marched quite well. 

The thermal structures of the two 1000-frequency models are generally 
Identical to within a few degrees, with a few differences as large as 10 # K. 

The differences are amiller than the computational and observational noise. 

The model computed with 2000 frequency points (with equal spsclng of 12. 5 cm * 
below 25000 cm’\ and 125 cm A above ) Is nearly Identical with the equally- 
spaced 1000-frequency model, llieae comparisons are not Illustrated since the 
models are Indistinguishable If drawn on the same acale as Figure 1 and 2. 

Figures 1 and 2 and the results of the previous paragraph strongly suggest 
the self-consistency of the OS treatment. Imagine the envelope of all models 
computed with different sets of N frequencies. As the number of frequency points 
N Is Increased, the resulting envelope of model atmospheres narrows nicely, until 
(near N • ’000 points) It defines a unique model structure. In addition, the 
envelopes determined for successively greater numbers of frequencies fall within those 




A 


specified by che coarser grids. Altogether these tests provide persuasive 
evidence of the convergence of the OS method In both of the ways outlined at 
the beginning of this section. 

It appears that 1000 fr ( uencles are sufficient to fix the structure 
of a model atmosphere within 30*K. This uncertainty Is very acceptable, 
even for fine analysis abundance determinations or synthetic spectra. In fact, 
the 500 -frequency models are sufficiently precise for moat purposes. This 
conclusion on the critical number of frequency points agrees well with that 
of Peytremann (1974). The only Important consideration In the selection of 
the frequency points la that they he distributed over the entire spectrusi. 

B. Comparison of Atomic and Molecular Line Blanketing 

The separate Influence of atomic and molecular lines on atmospheric 
structure is Illustrated In Plgure 2, which compares standard atmospheres 
(40CU/2.0/S /HOP, OS) computed with different amounts of hlanketlng. The curves 
correspond to: (1) continuum opacity only, (2) continuum plus atomic line 
opacity, or (3) continuum plus atomic line plus molecular line opacity. The 
effect of the atomic lines Is as expected: the nurface layers are cooled 

(f> 2400*10 In the blanketed model compared to the unblanketed model (which 
has 3250*K), and the deepest layers are warmed slightly. The additional 
blanketing of the molecular ltnja further cools the outer layers (to 2040*K 
and re ses the temperature of the deepest layers slightly more. For the 
standard model, most of the effect from the molecules Is due to CO. The strong 
surfce cooling effect and slight backwarmlng for CO has alread" been demonstrated 
through the use of straight mean opacities (Johnson 19731, nno Is confirmed In 


a more precise way here. 
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C. OS and Straight Mean Opacities 

Figure 3 Iso shows a curve for an atmosphere computed with straight 
nean opacities taken over Intervals of 100 cm *. Straight mean opacities, 
which are simply the average of the opacity across a chosen spectral interval, 
have the great advantage of simplicity and have been used by a number of 
investigators (for references, see Paper 1). The particular model shown In 
Figure 3 Includes a fraction of the MK atomic line blanketing (Mutschlecner 
and Kelli r 1972) and straight mean opacities of H 2 O, CO, and CN. While some 
of the differences between the straight mean and the OS models may be due to 
slightly different atomic and molecular data, most of the differences must be 
due to the more accurate treatment of the line opacities in the OS models. The 
principal error In the strlght mean opacity arises from the smearing of the opacity 
and the consequent filling of t».e opacity windows, through which more flux would 
otherwise flow. The straight mean therefore tends to overblarket a model, resulting 
In higher temperatures In the upper atmosphere and perhaps throughout the model 
(Carbon 1974), 

We caution, however, against uncritical Inferences about the straight 
mean method from this comparison. Like other such comparisons In the literature, 
we have here comparud a strlght mean model calculated with only 100 frequency 
points with an OS model computed with 1000. Clearly, the accuracy of the straight 
mean must depend on the size of the frequency Intervol over which the opacities 
are averaged. To our knowledge, the effect of the Interval size on the resultant 
model has never been carefully investigated, though some appreciation for the 
effect for the CN red system can be gained from the figures In Johnson, Marenln, 
and Price (1972). In the limit of an Infinitesimally small spectral Interval, 
the straight mean must he exact of course. The crucial question can then be 
stated as follows: What Is the critical size of the spectral interval over which 
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the straight mein la taken In order that the reaultant model approach the 
"correct" model within an arbitrary accuracy? 7l»ls would appear to he an 
Interesting area for further investigation. The question has recently been 
touched upon by Van Paradljs and Vardya (1975), who claim that even the 
opacity distribution function la a certain kind of straight mean. Perhaps it 
la more meaningful to think of the straight mea'. opacity as a one-picker ODF. A 
closer examination of the relationship betwen the straight mean, the OOP, 
and the OS method would appear fruitful. 

D. Comparison of the OS and OOF Model 

A crucial test of the ODS method would be the comparison of a model using 
the OS with one computed using the accurate opacity distribution function 
(OOF). This has not get ally been possible In the past because different 
Investigators have used different sets of llnec, differing dissociation energies, 
oscillator strengths, and co.apos i t Ions ; however, close agreement between the two 
methods has been achieved In a particular solar model (Peytremann 1974). Now 
Bell, et. al, (1975) have computed models using their giant line method of l'ne 
blanketing which Is an ODF method described In detail by Gustafsson, et, al. (1975). 
Using the precise data employed by these authors in computing their ODF models, 
we make such a comparison for a model (Bell 1974' with the parameters (4000/2.25/S). 
The results of this test are shown In Figure 4, As la apparent, the models 
produced by the two methods are In excellent agreement. This accord Is heartening 
on all counts since It verifies the validity both of the opacity sampling method 
and the giant line method. The alight differences which arise In the very 
outer layers of the atmosphere are not to he regarded as serious. Tn the outer- 
most atmospheric lnye a, the giant line method falls because the strongest lines 
are artificially truncated by the highest picket In the ODF, and the temperature 
plateau reached la slightly too hot. 
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The OS me t hod on the other hand, falls In the very outer layers because the 
energy balance there Is controlled by a very few strong lines and the stall stirs 
connected with these very few lines are poor. This Is of no Importance, however, 
because other more serious neglects (departures from 1.TF. and radiative equl l IbrlunO 
would likely make the results there Inapplicable to real stars anyway. The 
excellent agreement of both methods throughout the remainder of the stars 
provides substantial evidence that we have found the "correct" model. 

E. Effect o? Mne Wings on Models 

The above tests were performed with the assumption of the boppler 
absorption profile only for both atomic and molecular lines. When wings are 
added to the atomic lines as described by the usual Voigt function, the 

-U 

atmosphere Is very slightly warmed in the outer layers from < 10*K for t(R) > 10 
to as much as 100*K at t (Rosseland) ■ 10 * (where the atmosphere Is very 
uncertain, however). Deep backwarmlng Is negligible. In these calculations the 
total damping constant f * ffradiatlon) + f(Van der Manls). For simplicity we 
assume T (radiation)” 1y (classical), although the models are rather insensitive 
to the precise value used. With ffradiatlon' ■ 1.5 v (classical' the thermal 
structure la changed by less than 2*K anywhere. 

IV. CONCLUSIONS 

We have subjected an opacity sampling (OS) technique for calculating the 
opacities of large numbers of atomic and molecular lines to a number of tests, 
with very promising results for the method. Our conclusions follow. 

(1) For a giant stt r of solar composition anJ effective temperature of 
4000°K, we find that 1000 frequency points, spaced either randomly or uniformly, 
produce atmospheric models which are identical within 10“K at the worst point 
in the atmosphere. Changing the location of the frequency points or adding 
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Additional frequency point.- beyond this critical number produces no further 
changes In the atmosphere. In that *tm< It la "converged". A a«t of 300 
frequency points, picked randomly or equally spaced, produces an atmosphere which 
agrees with the above standard atmosphere within 30*K except In the extrsoe outer 
atmosphere, where the agreeae nt is within 100*K at the worst point. iospherlc 
models comput :<i with as few as 100 frequency points spread randomly across the 
entire spectrum agrees quite well (within 30*K) with the standard model throughout 
most of the photosphere hut are rather badly In error In the outer layers where 
the T(Kosseland) Is less than 0.001. 

(?) For a particular model with the parameters fi000/2.0/S! we find that 
the addition of atomic lines in the OS technique cools the boundary temperature 
from about 3250*K to ?400*K and warms the deeper layers of the photosphere 
slightly. The addition of molecular lines of CO, CN, CH, and C,, (of which CO 
Is tb- ly significant absorber for this model) cools the boundary to a 
tfSj-ernture of 2040*K and warms the deep photosphere even more. 

(3) In a precise comparison of atmospheric models produced with the 
opacity sampling fOS) method and the opacity distribution function (ODF) method 
for the parameters (4000/2.25/S) , we find the tenperature-pressure to be Identical 
to within 10*K at the worst point in the atmosphere. This test provides persuasive 
evidence that both ol these methods produce a "correct" atmosphere. 

(4) For models of giant stars, the lint! wings add no significant opacity 
to the atmosphere, except a slight warming In the outermost atmospheric layers. 

(5) The results of this paper strongly support the results of a previous 
paper (Paper I) regarding the advantages of the OS method. 

The authors have greatly benefited from helpful discussion* with C. neden, 
T. Fay, D. Carbon, G. Ryblckl, G. Collins, and R. Pall. Special thanks go to 
R. Bell ano R. Kuruce for their kindness in allowing use of their data tapes. The 
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Captions for Figures 

Temperature-pressure structure of a standard atmosphere 

(as la Fig. 1) converged with different sets ol 100 frequencies 

and one set cf 1000 frequencies (solid line). 

Temperature - pressure structure cf a standard atmosphere 
14000/2. 0/S) converged with different sets of S00 frequencies 
and one set of 1000 frequencies (solid line). 

Effect of atomic lines and of both atomic and molecular 
lines compand to continuous opacities alone on the structure 
of a standard model (as In Fig. 1). A model with straight mean 
opacities Is also given for comparison. 

Comparison of models (4000/2. 25/8) computed by the gl/ nt line 
(OOF) method and by the opacity sampling (OS) method for 
essentially the same set of atomic and molecular parameters. 
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